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O ■ Abstract 

^~>. The liquid and supercooled states of water show a series of anomalies whose nature is lively 
^ ■ 

, Qh | debated. A key role is attributed to the formation of structural aggregations induced by critical 

,__i ' phenomena occurring deep in the supercooled region, not experimentally accessible for bulk water. 

f"~» ■ This explain why, despite the numerous experimental investigations, the nature of water anomalies 

OO , and the hidden critical processes remain elusive. Here we present a time-resolved optical Kerr 

T^lj- ' effect investigation of the vibrational and relaxation processes in supercooled bulk water. The 

en ' experiment measures the water intermolecular vibrations and the structural relaxation process in 

ILJ . an extended temperature range, and with unpreceded data quality. A mode-coupling analysis of 



^ I the experimental data allows the characterization of the intermolecular vibrational modes and of 

H ■ 

their interplay with the structural relaxation process. 



C^ 



I. INTRODUCTION 

Water can remain liquid below its melting point and stays in its metastable phase known 
as supercooled water [l|,l2|. Several thermodynamic (e.g. isothermal compressibility, isobaric 
heat capacity and thermal expansion coefficient) and dynamic (e.g. viscosity, diffusivity, 
structural relaxation) observables of liquid water show an anomalous behaviour [5|]. In the 
supercooled state the values of these physical observables clearly diverge, with a singularity 
temperature Tg estimated at ^ 228 K. Despite the wide interest and the numerous research 



efforts on this subject, the origin of these critical phenomena remained unclear 



lively debated 



QI9 



and 



15|- 



It is widely accepted that the water anomalies are connected with the formation of var- 
ious structural aggregates whose features are modified by the temperature/pressure condi- 
5|. A crucial point is the physical origin of these structural aggregations. 



tions 



a 



In their pioneering work Poole, Sciortino, Essmann, and Stanley [3| postulated the ex- 
istence of a liquid-liquid critical point (LLCP), based on a simulation study and on the 
experimental observation of polyamorphism in glassy water (i.e. the existence of two forms 
of water glass: low-density and high-density amorphous ice). In the LLCP hypothesis a ffist 
order phase transition take places involving two supercooled water liquid forms: the low- 
density and the high-density phases. The low-density phase is characterized by tetrahedral 
intermolecular coordination, while the high-density one shows more packed structures with 
distorted networks. The existence of a LLCP is compatible with a series of computational 
and experimental results J5|. Recent new simulation studies con&med the LLCP hyphote- 
sis 1^, Il2l4l5|. Nevertheless, the matter is still under debate: in their recent simulation 



study, Limmer and Chandler interpreted the water anomalies as non-equilibrium phenom 



ena associated to the liquid- crystal transition 



Q 



Moreover the relevance of ice nuc 



eation 



phenomena and their relationship with the water anomalies has been reported |9l. Il0|. 

A fundamental question in order to explain the origin of water anomalies is if these two 
liquid forms persist at temperature higher than T^, where they must be considered as lo- 
cal and transient structures. The fluctuations between the two water species could explain 
the thermodynamic and dynamic anomalies. Unfortunately, the most evident critical phe- 
nomena are expected to occur in the supercooled phase at temperature/pressure values not 
directly accessible in bulk water samples, and this prevents a direct experimental solution 



of the problem. In fact, in several experimental and simulation works the water singularity 
temperature Tg at atmospheric pressure has been extrapolated to be around 223-228 K, while 
homogeneous nucleation fixes the lowest reachable temperature at about 231 K for liquid 
water. Actually, in a macroscopic bulk sample the crystallization phenomena and several 
experimental difficulties limit the minimum temperature around 243 K. Nano-confinement 
has the advantage of preventing water freezing and of extending the experimental investiga- 
tion of the supercooled phase down to very low temperatures; it is however still unclear to 
what extent the interactions with the pore surfaces modify the water properties 16|-[l9|. The 
hydrogen bond network largely determines the intra- and inter-molecular vibrational spec- 
trum of supercooled water; thus the experimental investigation of the spectral features gives 
precious information on the local intermolecular structures and on their variation with the 
temperature/pressure conditions. The investigation of high frequency spectrum, 1000-4000 
cm~^, gives access to the intra-molecular vibrational dynamics. In a recent paper Mallamace 



et al. (20[ investigated the OH vibrational spectra of confined supercooled water by means 
of Fourier transform infra-red spectroscopy (FTIR). According to their experimental results 
a new vibrational component at 3120 cm~^ appears at low temperature, whose intensity 
increases sharply in the deep supercooled phase. The authors attribute this component to 
the low-density water liquid form. The low frequency spectrum of water, u < 1000 cm~^, 
refiects the inter-molecular dynamics and so it is particularly sensitive to the local molecu- 
lar structures. This frequency range of the water spectrum has been investigated in several 



light scattering studies, see for example 2ll-l27l|. Two main broad peaks are observed at 



50 and 200 cm~^, generally attributed to the hydrogen bond "stretching" and "bending" 
vibrations. The spectra show also a very low frequency wing, z/ < 20 cm~^, that has been 
attributed to "relaxation" processes. The light scattering data were measured at relatively 
high temperatures and have rather poor signal-to-noise ratio, due to the very weak water 
signal: these drawbacks did not allow identifying all the possible vibrational modes possi- 
bly contributing to the water spectrum in this region. Moreover, the reported attempts of 
fitting the light scattering data have been based on a simple superposition of independent 
Gaussian and Lorentzian functions, neglecting any possible coupling between the vibrational 
dynamics and the structural relaxation processes. 

An accurate analysis of the low frequency spectral features of supercooled water, includ- 
ing a quantitative comparison to the predictions of recent models for water dynamics is 



still missing; it actually appears as a required step towards the full understanding of the 
structural and dynamical properties of water in its low temperature metastable state. 

In recent years, time-resolved non-linear spectroscopy demonstrated to be a very use- 
ful tool for the investigation of complex liquid dynamics 28|. Several studies of this type 
have been performed on liquid water; only in few cases they have been extended to its su- 
percooled phase 29l-l32|. Time-resolved optical Kerr effect (OKE) investigations of liquid 



water 



33 



show the signature of fast vibrational dynamics at short times, followed 
by a slower monotonic relaxation. The initial oscillatory component provides information 
about intermolecular hydrogen-bond dynamics, and can be directly compared with the low- 
frequency Raman spectrum 33|. Actually, the OKE experiments on supercooled water 29 1 
have been focused so far on the investigation of the slower decay, which is interpreted as a 
structural relaxation phenomenon. No detailed OKE investigation of the vibrational dynam- 
ics in supercooled water has been carried out, due to the intrinsic experimental difficulties. 



so that a detai 
mental set-up ■ 



ed picture of this dynamics is still missing 37|. Using an improved experi 



3l|, we were able to measure, for the first time according to our knowledge, 



the fast vibrational dynamics and the slow structural relaxation in a single data set char- 
acterized by low noise and wide dynamic range. Moreover, we reached temperatures lower 
than in previous experiments, approaching the homogeneous nucleation limit. We succeeded 
in measuring the OKE signal relaxation over a large time interval of about 25 ps with a high 
time resolution of 20 fs. We can thus provide the unambiguous measurement of the entire 
correlation function in supercooled water, enabling a very stringent comparison with the 
available theoretical models. 

II. EXPERIMENT AND METHODS 



Time resolved optical Kerr effect (OKE) spectroscopy is a non-linear technique exploiting 
the transient birefringence induced by the electric field of radiation in an optically trans- 
parent and isotropic medium. The method has been described in many papers and review 



articles 



28 
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-40|. In short, the oscillating electric field of a polarized short pulse of 
monochromatic light (the pump) induces optical birefringence in the sample. In addition to 
the instantaneous electronic contribution (which does not carry any dynamic information), 
the induced birefringence reflects the relaxation and vibrational response of the molecules in 



the sample. A second pulse (the probe), spatially superimposed within the sample with the 
pump pulse and characterized by a different polarization, probes the induced birefringence. 
The intensity of probe light leaking through a polarizer placed after the sample measures the 
material response as a function of the time delay between pump and probe, the optical set- 
up is reported in the Supplementary Information. If femtosecond light pulses are employed, 
the method allows measuring the relaxation processes and the low frequency vibrational 
correlation functions of the system. In this sense it is the time-domain counterpart of de- 



polarized light scattering technique |40|. The wide accessible time window, extending from 
tens of femtoseconds to hundreds of nanoseconds, makes OKE a very flexible technique, able 
of revealing very different dynamic regimes, from simple liquids to supercooled liquids and 



glass formers 29|, l4l|-|43|. In the case of optical heterodyne detection, the delay time de- 
pendent signal is directly proportional to the material response (relaxation) function, R{t), 
convoluted with instrumental function, G{t). The response R{t) is directly connected to 
;he time derivative of the time-dependent correlation function of the dielectric susceptibility 



m 



m 



m ex |(x(t)x(0)) (1) 

The Fourier transform of R{t) correspond to the frequency- dependent response, the observ- 



able measured in depolarized light scattering (DLS) experiments 28|,|40|]. We stress here that 
measuring the correct instrumental function G{t) is a crucial step for achieving a reliable 
and accurate fitting of the experimental OKE signal. A detailed description of the method 
adopted is given in the supplementary information. In short, the critical point in the entire 
procedure is that any minor, apparently negligible adjustment of the optical set-up when 
passing from recording the instrumental function to measuring the water signal prevents the 
extraction of the real instrumental function. The resulting inaccurate deconvolution of the 
water signal degrades the quality of the overall fitting. In our experiments G{t) was ob- 
tained by measuring the OKE signal of a CaF2 plate placed in a purposely designed holder 
which allowed switching to the water sample without even touching the rest of optical set- 
up. Supercooling bulk water can be achieved only with particularly pure sample and using 
specific attention for the cooling procedure. We performed the measurements on a sealed 
vial of cylindrical shape, prepared for pharmaceutical purposes by the Angelini company, 
which allowed us to reach the temperature of 247 K. The sample temperature is controlled 
by means of home-made cryostat with a stability of 0.1 K, using a Peltier cooler. The laser 
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FIG. 1. HD-OKE signal as function of temperature in log-log plot. At short delay times (<1 ps), the 
HD-OKE signal shows oscillations due to vibrational dynamics, that become more defined as the 
temperature is lowered below the melting point. At longer delay times, the signal is characterised by 
a monotonous decay, that has a strong temperature dependence and a clear non-exponential nature. 
For a readable display, only 6 temperatures (247-258-273-293-313-353 K) out of the 11 measured 
(the previous ones plus 253-263-268-278-333 K) are reported; the intensities are normalized at 
0.7 ps. The upper and lower insets show, in a log-linear plot, the 247 K HD-OKE data at short 
and long delay times, respectively. 



system is a sel: 



-mode-locked Ti:sapphire laser producing pulse of 20 fs duration with an 



energy of 3 nJ 31|. The optical set-up and other details about the experimental apparatus 



are reported in the supplementary information. 

III. RESULTS AND ANALYSIS 

The heterodyne-detected OKE (HD-OKE) data are reported in Figured) the data clearly 
show the complex dynamic features presented by liquid and supercooled water. The vi- 
brational dynamics, extending up to 1 psec, displays an articulate pattern that becomes 
progressively more defined going deep in the supercooled phase, and that at longer times 
merges into the monotonic decay due to structural relaxation. The time-resolved HD-OKE 
experiment measures the time-dependent correlation function of the collective water polar- 



6 



izability, see Equation ([T]). The time evolution of this function reflects the water dynamics 
thorough the electronic polarizability. Numerical investigations 37| show that the OKE 
signal of water, in all time scales, is the results of a complex interplay between intrinsic 
molecular terms and interaction-induced contributions. In particular, the fast oscillatory 
dynamics consists of a combination of collective hindered roto-translational motions. 

The relaxation features presented by supercooled water are very similar to that of the 



glass-former liquids 43|, and the slow dynamics is found in agreement with the main mode- 



coupling theory predictions 29|. In Figj2] we report the master plot analysis of the slow 
component of the HD-OKE signal decay. The new data confirm, over an extended temper- 
ature range, the power law temperature dependence of the relaxation time. The analysis of 
the whole time- dependent correlation function measured by the optical experiments requires 
the utilization of a model that enables an operative parametrization of the vibrational modes 
present in the wate r dy namics. This turns out to be a complex task, and different attempts 



have been pursued |34|-l3q| in that direction. None of these water studies was extended into 
the supercooled state. Moreover, the data analysis and fitting models adopted so far did 
not include any coupling between the fast vibrational dynamics and the slow structural re- 
laxation. We tested some of these fitting procedures on our HD-OKE data and found that 
they fail to reproduce correctly the HD-OKE data in the supercooled temperature range. 



In order to analyze the vibrational dynamics and the structural relaxation inc 
coupling processes, we utilized the schematic solution of mode-coupling theory 



uding their 



44| 



The formalism of mode-couplin g th eory (MC) stems from a generalization of the Mori 
and Zwanzig hydrodynamic theory 4^. In MC theory the retardation effects are taken into 
account by memory functions M(t), which play the role of the hydrodynamic friction coeffi- 
cients. Of course the task of any microscopic theory is that of providing the form of M(t). In 
its simplest form, developed for simple atomic liquids, only density fluctuations contribute 
to memory; in order to overcome the limitations of this approach and to extend the MC 
theory description to the dynamics of more complex systems, several versions of the theory 
have been developed. For molecular liquids in fact orientational and intra-molecular degrees 
of freedom contribute to the overall dynamics. In our analysis we adopt the "schematic" 
formulation of the MC theory (SMC). Details of the theory are given in the supplementary. 
Briefly, in this approach the density-density time correlator, $,m oc (p(t)p(O)), termed the 
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FIG. 2. The slow dynamics of supercooled water can be analyzed using a model-independent pro- 
cedure 291] , consisting in rescaling the time/amplitude scales of the raw data measured at different 
temperatures in a way that the decay curves overlap. If all the data collapse in a "master-plot", it 
is possible to identify a function, the "master function" , that defines the relaxation independently 



of the temperature 



44 ] ■ In the upper panel we report the master-plot of the HD-OKE data and the 



master function: ■^exp{—x^) with a stretching parameter /3 = 0.6. The relaxation times obtained 
with the rescaling procedure are shown in the lower panel of the figure. The new data (red circles) 



291 ] ■ The temperature dependence of these 



extend and confirm the previous data (empty circles) 

data are in agreement with the power law: r oc (T — Tc)~'^, with a critical temperature Tc — 227 K 

and a critical exponent 7 ~ 1.7. We report the relaxation times and the power law in a linearized 



plot, in order to visualize immediately the scaling proprieties and the critical temperature. 



master correlator, is described by the differential equation: 

$m(t) + Vm^mit) + ^^^^(t) + f Kit- t')K.{t')dt' = 0, (2) 

where Qm is the mode frequency, rjm is the shear viscosity and the memory function K{t) is 
expressed as a series expansion of the master correlator itself. The role of other observables, 
Ai, linked to the density correlator dynamics (e.g. the molecular orientation, the dynamics 
of a target molecule, etc.) is described by a set of differential equations [39] of the same form 
as equation (|2]), written in terms of slave correlators $i(t) oc (Aj(t)Aj(O)) and of memory 
functions of the form mi{t) = v^^rn{t)^i{t), which implies the coupling of the master and 
slave correlators. The key point in the SMC analysis of the OKE experimental data is that 
the time correlation function of the sample polarizability can be expressed as the sum of the 
time derivatives of the slave correlators: 

i?^^-^^(t)oc|$:$.(t). (3) 

i 

We underline here that the individual slave correlators describe "average" collective modes, 
so that they cannot be associated to specific physical variables. In summary, the SMC 
approach represents a robust physical model able of describing a complex dynamics, with 
no need of unjustified assumptions concerning the decoupling of slow and fast dynamics. 

We solved the SMC equations numerically, taking the frequencies, viscosities and coupling 
coefficients as parameters to be adjusted in order to reproduce the HD-OKE response by a 
fitting procedure. 

The details on the SMC equations and fitting procedures are reported in the supplemen- 
tary information. In Figure |3] we show some of the measured HD-OKE signals with the fits 
obtained using the SMC model. The model is able of reproducing correctly the experimental 
data over the whole time range at all temperatures. Most remarkably, and differently from 
other fitting models widely used, this result is achieved without imposing any non-physical 
decoupling of vibrational and relaxation dynamics. An important result of our work is that 
the fitting at temperatures < 313 K requires three slave correlators, whereas for higher 
temperatures only two correlators are needed. 

The time domain data, collected in the time-resolved HD-OKE experiments, are trans- 
formed to yield the corresponding spectra in the frequency domain. In fact, the Fourier 
transform of the HD-OKE response function, R{u}), is directly connected with the spec- 
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FIG. 3. The complex vibrational dynamics taking place in the sub-picosecond time scale is correctly 
reproduced by the numerical solution of the SMC model. The SMC model enables a very good 
fit of the slow relaxation decays, corresponding to the stretched exponential function. Also at 
intermediate delay times, where the vibrational dynamics merge into the relaxation processes, the 
decay curve is correctly described by the SMC equations. We report the SMC fits (red line) of the 
HD-OKE data (circle) for five temperatures. In the left panel the short time window is shown in 
a log- linear plots, in the right panel the whole time window is shown in a log-log plots. 



trum Sdls{^) measured in the depolarized light scattering (DLS) experiments: the imagi- 
nary part of R{uj) is proportional to Scisi^) corrected for the Bose factor, i.e Im[R{u)] oc 



Sdls{^)/[''t-{^) + 1] (28|,|40||. Our results then can be directly compared to the DLS, Rayleigh 



and Raman, found in the literature. The imaginary part of the Fourier transform of the 
HD-OKE data and the calculated response are reported in FigJH The experimental spectra 
clearly show the characteristic vibrational features of liquid water: the two broad bands 
around 50 and 200 cm~^, and the low frequency wing. This wing corresponds to the slow 
decay of the OKE signal measured in the time domain; as discussed in our previous work 29| . 
it is attributed to structural relaxation. 



10 



The other two bands are generally attributed to translational intermolecular modes (i.e. 
the oscillating dynamics of water molecules around their instantaneous center of mass) hin- 



dered by the first neighbour molecular cage 



37 



45 



46| . In the liquid phase they are largely 



affected by the hydrogen bonds. A comparison with the calculated vibrational density-of- 



states reported in ref. 45|] suggests specific assignment of the 50 and the 200 cm ^ bands. 



The first one is dominated by the transverse translational motion, while the second has essen- 



tially longitudinal character 45|. The physical interpretation is that the two intermolecular 



vibrations are mainly involve "bending" and "stretching" distortion of the hydrogen bonds, 



respectively 23| . This picture even if intuitively useful, is clearly oversimplified: it doesn't 
take into account the coupling of the two vibrations imposed by the complex hydrogen bond 
network present in the liquid water. The presence of two different local structures, the 
one corresponding to highly tetrahedral intermolecular structures, and that consisting of 
strongly distorted network, as well as their continuous reorganization induced by the struc- 
tural relaxation processes, are expected to strongly affect the vibrational dynamics. The 
SMC fit enables to disentangle the vibrational modes contributing to this dynamics and to 
point out their coupling to the structural relaxation processes. The bending mode at about 
50 cm"^ is described by a single correlator (magenta line in FigJl]) that turns out to be 
decoupled from the structural relaxation (the peak close to frequency). In contrast, the 
stretching mode, around 200 cm~^, needs two correlators to be properly reproduced. One 
correlator, peaked at about 180 cm~^ (blue line), shows an almost temperature independent 
peak amplitude and a large spectral component extending to very low frequency, indicative 
of strong coupling to the structural dynamics. The other correlator, peaked at about 220 
cm~^ (green line), appears to be substantially uncoupled to the structural dynamics; its 
amplitude increases with decreasing temperature. Our HD-OKE measurement represents 
the first clear experimental determination of this feature, even if it was inferred in a previous 



Raman study |21 |. 



A comparison between these results and the recent debate on the nature of water anoma- 
lies is due. Our data and analysis show the presence of two distinct modes responsible for the 
part of the water spectrum normally attributed to the inter-molecular stretching vibrations. 
Raman spectra of low and high-density amorphous ice (LDA and HDA) in this frequency 



region differ from that of liquid water 47|]. LDA has a broad peak centred at higher fre- 



quency, about 215 cm ^, HDA shows an even broader peak centered at about 185 cm ^ 

11 
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FIG. 4. The Fourier transform of the time-dependent HD-OKE response, -R(w), can be directly 
compared to the frequency-dependent spectrum measured in the DLS experiments, SoLsi^)- The 
imaginary part of the Fourier transform of the HD-OKE data (o), after that the contribution of the 
instrumental function has been removed by de-convolution (see supplementary information) and of 
SMC fits (red line) at 247-258-273-313 K are reported in the frequency scale, v = uj/2Tr. The SMC 
model allows to disentangle the dynamical features that generate the response function used in the 
fit. The contributes of the three slave correlators are reported in the figure as magenta-blu-green 
lines. The contribution of the very high frequency components {u > 400 cm~^) of the response 
is included in the fit by means of a simple damped harmonic oscillator. This does not affect the 
response function in the lower frequency range, where the dynamics relevant to our investigation 
takes place. 
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(see for example fig.l in ref. [47|) whereas the crystalhne ice phases typically have a narrow 
peak at 230-240 cm"^. This suggests that the modes observed in our experimental investiga- 
tions can be associated with two fluctuating water species with different local structures; a 
low-density form (LD) characterized by tetrahedral network and four-coordinated molecules 
similar to the LDA structure, and a high-density (HD) form characterized by closely packed 
aggregates with lower coordination and high network distortions. The mode at 220 cm~^ 
can be attributed to the LD form and the mode at 180 cm~^ to the HD one. The present 
data strongly support the existence of the LD and HD water species. Nevertheless, as the 
observed water djTiamic is measured at temperature/pressure relatively far from the sup- 
posed LLCP, some ambiguity remains about the nature of the critical phenomena taking 
place into the supercooled phase. The spectroscopic features attributed to the LD and HD 
forms could be attributed either to the coexistence of two extended liquid phases leading to 
the liquid-liquid phase transition (l|, y, |5|, or to the continuous transition from a normal, 



Iq 



unstructured liquid to a tetrahedrally structured liquid J9|. 



Our new experimental data and analysis clearly show the presence of the high-frequency 
satellite mode contributing to the intermolecular stretching band of water, whose amplitude 
increases sharply when the temperature goes deep in the supercooled phase. The nature of 
this mode, peaked at 220 cm~^, appears clearly different from that of the nearby mode at 
180 cm"^. In fact, the former is strongly temperature dependent and only weakly coupled to 
the structural relaxation processes, while the latter is almost temperature independent and 
strongly coupled to structural rearrangements. The surprising similarities with the spectra 
of LDA and HDA glass forms with that of the two high frequency distinct modes, reported 
in this work, suggest attributing them to the LD and HD water structures, respectively. The 
LD water forms are expected to consist of aggregates of a small number of molecules; being 
mostly localized in space (at least in the temperature range investigated), their dynamics 
show negligible coupling with the structural relaxation. The number/lifetime of the LD 
species must increase with decreasing temperature in order to account for the amplitude 
increase of the corresponding correlator revealed by our analysis. On the other hand, the 
vibrational dynamics of the fraction of HD water appears dominated by the bare caging 
effect that merges into the structural relaxation, a behaviour typical of fragile glass-former 
liquids. 

13 
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Appendix A: Supplemenraty Information 
1. Experimental features 

Supercooling bulk water is not an easy task because of crystallization. Temperatures as 
low as 20 degrees below melting can be achieved only with particularly pure sample and 
using specific attention for the cooling procedure. In our work we performed the measure- 
ments on bi-distilled water prepared by the Angelini company in a sealed vial of cylindrical 
shape, which enables us to reach a temperature as low as 247 K. The sample tempera- 
ture is controlled with a stability of 0.1 K by means of a home-made cryostat. The vial is 
inserted into a parallelepiped-shaped aluminium holder, with a central cylindrical housing 
whose diameter is close to that of the vial. Two fused silica windows on opposite sides of 
the aluminium holder allow the beams to cross the sample. A thin film of glycerol between 
the vial and the housing assures an efficient heat transfer. The holder is fixed to the cold 
plate of a Peltier cooler and the temperature is controlled by a platinum thermoresistance 
in thermal contact with the holder itself. 

A self-mode-locked Ti:sapphire laser (by Femtolasers, mod. Fusion: 25 fs pulse dura- 
tion, 3 nJ pulse energy, 77.47 MHz repetition rate, 800 nm wavelength) was used for the 
experiments. Water is a liquid that has a very low OKE signal due to its nearly isotropic 
molecular polarizability. Our set-up enables measurement of the fast vibrational dynamics 
and the slow structural relaxation in a single experiment, with high signal-to-noise ratio. 
The delay time was scanned by means of the continuous motion of a translation stage, 



whose absolute position was controlled by a linear encoder 48|. This method reduces the 



acquisition time substantially, while improving the signal statistics. We adopted a peculiar 
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FIG. 5. Heterodyne detection scheme of OKE esperiment. 



configuration for the heterodyne detection SjJ, |49| , with circularly polarized probe beam and 
differential acquisition of two opposite-phase signals on a balanced double photodiode. With 
reference to FigJSl a quarter wave-plate placed between the two polarizers produces a circu- 
larly polarized probe field. Two signals, corresponding to the horizontal and vertical linear 
polarizations generated by polarizer P2 and having a opposite sign, are sent to the balanced 
photodiode detector. Thus the measured OKE signal is automatically heterodyned and free 
from possible spurious signals. A further improvement of the data quality was obtained by 
subtracting the two HD-OKE measurements, corresponding at left and right circular polar- 
izations of the probe field. This procedure removes from the signal the spurious dichroic 
contributions coming from possible misalignment of the wave-plate. 



2. HD-OKE signal 



The signal measured in an HD-OKE experiment is 50|, l51| : 

S{t) = f dt\Eprit - t)\^ fdt'R{t-t')\E,,{t 



'm2 



(Al) 



where Ep^ and E^x are the probing and exciting laser fields, respectively; R{t) is the material 
response function. Since we are performing a non-resonant OKE experiment, the Born- 



Oppenheimer approximation applies and the response function can be cast in the form 51 
53|: 

R{t) = ^6it) + Rnit) (A2) 

with 7 representing the instantaneous electronic response and -R„(t) the nuclear response. 
The latter can be rewritten as 

Rnit) = -lll^xx (A3) 
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with $^^ = {Xxy(t)Xxy{0)) where Xxy(t) is the off-diagonal element of the collective electronic 
polarizability tensor. The Fourier transform of time-dependent response, R{u}), can be 
directly connected with the signal Scisi^) measured in the depolarized light scattering 
(DLS) experiments. In fact, the imaginary part Im[R{uj)] is proportional to the DLS signal 



measured in forward configuration, i.e Im[R{uj)] oc Sdls{^)/[''^{^) + 1] [51 1. 

In order to fit the OKE data we simulated the measured signal using the following ex- 
pression |38|, |5J] : 

S{t) = f [-fSit - t') + Rn{t - t')] G{t')dt' (A4) 

where G{t) is the instrumental function accounting for the finite time resolution determined 
by the time width of the pump and probe pulses. Under specific conditions [38|, G{t) 
corresponds to the temporal autocorrelation of the laser pulse intensity. The experimental 
investigation of the instrumental function is the argument of the following paragraph. 

3. The instrumental function measurement 

An important experimental issue concerns the measurement of the real instrumental 
function G'(t), which is far from trivial. By definition, the instrumental function in a HD- 
OKE experiment is the signal measured in a system with response function R{t) = S{t), 
as directly inferred from eq JA4[ In practice, this condition can be fulfilled by choosing a 
sample with negligible nuclear response. In most of the papers published on this subject in 
recent years a fused silica plate has been used to this purpose, as its nuclear contribution is 
considered negligible. 

In reality, silica presents a nuclear contribution, known to the glass physics community as 
"Boson Peak" . Although this contribution is very weak, it prevents from measuring the real 
instrumental function. This is particularly relevant for water measurements, especially when 
investigating the fast dynamics. Alternatively, the instrumental function has been obtained 
by measuring the second harmonic autocorrelation of the laser pulse, or more simply by 
fitting the electronic contribution to the OKE signal with an analytic peak function (like 
Gaussian, hyperbolic secant, etc.). 

In principle the second harmonic autocorrelation should be the best way to measure 
the instrumental function. It directly gives the temporal autocorrelation of the laser pulse 
intensity and, consequently, the temporal profile of G(t). The procedure, however, implicates 
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some changes and adjustments in the OKE set-up: for instance, i) the insertion of a second 
harmonic crystal in place of the sample, with the consequent change of the pulse compression 
status, ii) the re-alignment of the beams and iii) the translation stage re-positioning in order 
to achieve the optimal temporal and spatial superposition of the pump and probe. All 
of these adjustment potentially modify the response function. Actually, any even small 
changes of the experimental conditions (e.g. a different beam alignment or a modified angle 
of superposition of the beams inside the sample) affect the instrumental function. Therefore, 
any method for obtaining the instrumental function that involves changes of the set-up is 
unsafe and can yield an incorrect time profile. 

On the other hand, the accurate knowledge of the real instrumental function is mandatory 
to extract the true OKE response of water in the fast dynamics range. For this reasons we 
decided to extract the instrumental function by performing HD-OKE measurements in a 
reference sample. This procedure, differently from the second-harmonic approach, allows 
us to accurately preserve the experimental conditions used during the water measurements. 
As the reference sample we chose a plate of calcium fluoride {CaF2). This a cubic ionic 
crystal which is optically isotropic and, in the frequency range probed in our experiment, 
has only one Raman active band at 322 cm~^, due to the optical phonon T2g symmetry. So 
the nuclear response of the reference sample is a simple and well known functional form. 
The calcium fluoride plate was dipped inside a water vial, identical to that used for water 
measurements, supported by the same sample holder. 

The measurements of the instrumental function were done by just replacing the water vial 
with the water-CaF2 vial, leaving the rest of the set-up unchanged. In the measurement we 
took care that the faces of the CaF2 plate were perpendicular to the line bisecting the angle 
formed by pump and probe beams. The thickness of the CaF2 plate, 3 mm, was enough to 
fully contain the probe and pump overlapping area, avoiding any spurious signal contribution 
by the outer water. We took a reference measurement for each set of water data. Maximum 
attention was taken to preserve the experimental conditions between the water and the 
reference measurement. Figure reports a typical HD-OKE signal obtained in the CaF2 
reference sample: the signal shows a first peak, characterized by a very fast rise and fall, due 
to the electronic response, and a second oscillating contribution due to the nuclear response. 
This signal is the convolution of the instrumental function with the OKE response, see 
eq JA4|, where the nuclear part is a single harmonic damped oscillator (DHO). The simplicity 
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FIG. 6. HD-OKE signal of CaF2 (circles). The continuous line is the instrumental function 
obtained with the fitting procedure described in the text. The nuclear response, Rn, is modelled 
as the time-derivative of a single damped oscillator and the instrumental function, G, is simulated 
as the sum of Gaussian, Lorentzian, and hyperbolic secant functions. The comparison with the fit 
of the instantaneous electronic part performed using just a single hyperbolic secant (dashed line) 
is reported too. 



of this nuclear response function enables a reliable extraction of the instrumental function, 
G{t), by means of an iterative fitting procedure (i.e. least square fitting of the HD-OKE 
data of the reference sample with the simulated signal according to eq JA4p . We found that 
the instrumental function can not be reproduced by a single hyperbolic secant function, as 
often reported in the literature; a good fit requires the sum of several analytical functions, 
typically a combination of Gaussian, Lorentzian and/or hyperbolic secant functions. In fig. [6] 
we report the instrumental function obtained by the iterative fitting procedure (continuous 
line) and the one obtained by fitting the electronic peak with a simple hyperbolic secant 
function (dashed line). Apparently these two instrumental functions are very similar but 
the small differences cannot be neglected for an accurate investigation of water fast OKE 
response. This is mostly evident when the data are Fourier transformed to the frequency 
domain. In fact obtaining the imaginary part of the frequency dependent response from 
the time domain data involves the division by the Fourier transform of the instrumental 
function, i.e. Im[R{uj)] oc Im{FT [S{t)] /FT [G{t)]}. Fig. [7] shows the comparison between 
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FIG. 7. Comparison between the calcium fluoride frequency dependent response calculated using 
the instrumental function obtained by the iterative fitting procedure (continuous line) and the one 
obtained by fitting the electronic peak with a hyperbolic secant function (dashed line). Only the 
knowledge of the real instrumental function allows measuring the correct response. 



the response calculated using the instrumental function obtained by the iterative fitting 
procedure (continuous line) and the one obtained by simply fitting the electronic peak with 
a hyperbolic secant function (dashed line). We clearly see that only an accurate measurement 
of the real time domain instrumental function yields the correct response. 



To our knowledge, this aspect of the OKE data analysis has never been addressed in 
such a quantitative way, and the present level of accuracy in analyzing the role the of the 
instrumental function in HD-OKE experiments is unprecedented in recent literature. In 
particular, we want to stress that the knowledge of the correct instrumental function is 
fundamental in order to extract the true OKE response, especially for the short delay time 
range. This is particularly critical for weak signals characterized by complex relaxation 
dynamics, as in the case of water. 
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4. Analysis by Schematic Mode- Coupling Model 



The schematic mode-couphng (SMC) equations can be written as follows |44j. ISSHSTj : 

^mit) + Vm^M + ^l^mit) + / K {t - t')^Ui')dt' = (A5) 

where the memory function is 

K{t)=V,<^m{t)+V2^l{t) (A6) 

The SMC model defines the memory by a series expansion on the master correlator itself 
$m and provides in this way a closed form for the integro-differential equation IA5I The 
SMC model identifies the master correlator as the density correlator $„ oc (p(t)p(O)). The 
dynamics of any other observable, Ai, linked to the density correlator dynamics (e.g. the 
molecular orientation, the target molecule dynamics, etc.) can be described by a similar 



differential equation [58l |: 

^,{t)+7]A{t) + n^,Ht)+ fm,{t~t')^,{t')dt' = (AT) 

where 

mi{t) = vt^mmiit) (A8) 

and $i(t) oc (Aj(t)Aj(O)) is the slave correlator. The coupling between the slave and master 
dynamics is assured by the product between the slave and master correlators in the memory 
kernel rrii. Finally, the experimental observable can be expressed as sum of these slave 
correlators: 

R'''^{t)^§-^J2Mt). (A9) 

i 

According to eq JA9l in our analysis of the HD-OKE data we are decomposing the elec- 
tronic polarizability correlator, ^^x' ^'^^'^ ^^^ ^^"^ '^^ ^ii't) correlators. It is implicit that 
the individual slave correlators cannot be directly associated to specific physical variables. 
Each of these correlators in fact describes an "average collective mode" whose dynamics 
is described by the SMC equations. The vibrational and relaxation features and coupling 
processes are included in the SMC equations by definition. In this respect the present SMC 
equations IA5tlA9l represent a robust physical model able of describing a complex dynam- 
ics including vibrational and structural relaxations, implicitly accounting for their mutual 
coupling. Differently from other approaches, the method does not require any decoupling 
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or dynamic separation between the fast (vibrational) and the slow (relaxation) processes. 
Thus, once the temporal relaxation function of the master correlator is known, it can be 
used to calculate the relaxation of the slave correlators and, eventually, that of the OKE 
signal. 

We found that the OKE data can be described in the entire temperature range by em- 
ploying three slave correlators at most. In particular, the weight of the highest frequency 
contribution decreases monotonically with increasing temperature, becoming negligible at 
the two highest temperatures, where only two slave correlators are sufficient to reproduce the 
data. The parameters of the model are sixteen: the master equation parameters 77^, ^m, ^i 
and V2, the slave equations parameters ?7j, fij, v^ with i = 1, 2, 3 and finally the three ampli- 
tudes Oi, 02, and 03. The integro-differential equations IA5I and IA7I were solved numerically 
with a step by step second order Runge-Kutta algorithm. Preliminary series of fits provided 



a qualitative estimate of the temperature be 
similar analyses reported in the literature [54 



lavior of t 



m 
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le parameters. In agreement with 



62|, in the subsequent step we forced 



some of the parameter either to assume fixed values or to follow simple temperature trends. 
In particular, we found that the frequency Qm and the vertex vi of the master oscillator are 
almost constant in the whole temperature range, and we locked their values to 66 cm^^ and 
0.33, respectively. The second vertex V2, instead, was found to increase almost linearly with 
decreasing temperature, and it was forced to obey the linear dependence V2 = Q — 0.014T. 
The friction parameter rjm and the remaining parameters of the slave oscillators were left 
free of varying. In figure E] we show the temperature dependence of the slave frequencies 
Qi and the friction parameters 77^ and rji, and of the three vertices v^. The temperature 
behavior of the fitting parameters is consistent with the results obtained by the previous 



studies in glass-former liquids 54N56 



59-62 
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FIG. 8. Temperature behavior of the frequencies of the slave oscillators (circles, squares, and 
triangles refers to z = 1, 2, and 3, respectively) is reported in panel a). In panel b) we show the 
temperature behavior of the master and slave friction parameters (diamond refers to the master 
oscillator, circles, squares, and triangles refers to i = 1, 2, and 3, respectively). Temperature 
behavior of the vertices vf (circles, squares, and triangles refers to z = 1, 2, and 3, respectively) 
are reported in panel c). 
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